NON-LINEAR EVOLUTION AND PARTON DISTRIBUTIONS 

AT LOW X. 



M. LUBLINSKY 

Department of Physics, Technion - Israeli Institute of Technology, Haifa 32000, 

ISRAEL 

E-mail: malQtechunix. technion. ac.il 

We suggest a new procedure for extrapolating the parton distributions from HERA 
to much higher energies. The procedure suggested consists of two steps. First, we 
solve the non-linear evolution equation. Second, we introduce a correcting func- 
tion for which we write a DGLAP-type linear evolution equation. The nonlinear 
equation is solved numerically and estimates for the saturation scale, as well as for 
the gluon density at THERA and LHC energies are made. 



1 Introduction 

In the present talk, we present a new approach for extrapolation of the parton 
densities to low xo. By low x we mean values up to 10^'', which correspond 
to the energy range potentially covered by the THERA and LHC. 

A conventional approach to parton density evolution is based on the linear 
DGLAP equation which describes the gluon radiation leading to an increase 
in the number of partons. However, at low x the parton cascade become dense 
and the recombination processes start to be important. We are convinced that 
such system cannot be described by a linear evolution any more. 

Our motivation is based on two main problems of the DGLAP evolution. 
First, it predicts a very steep growth of parton distributions at low x violating 
the unitarity constraints. The second problem is common for all perturbative 
series which are asymptotic. In application to DIS processes it leads to twist 
OPE break down at low x, when the high twists become of the same order as 
the leading one. We believe that a non-linear evolution is a solution to both 
problems! The non-linear evolution accounts for the saturation effects due to 
high parton density and it sums high twist contributions. 

2 New approach for extrapolating the parton distributions 

In the colour dipole picture the evolution is applied to the imaginary part of 
the dipole elastic scattering amplitude N{r±,x\ b) for the dipole of the size rx 
elastically scattered at the impact parameter b. The amplitude N is the major 
unknown to be determined. Our approach to the problem is based on two steps 
in which we obtain as a sum of two terms: N — N + AN. As the first step. 
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Figure 1. Solutions of the Eq. (|l|) as a function of distance. 



the function N is found as a solution of the non- linear evolution equation (|T|). 
This equation is valid in the leading ln(x) approximation of QCD. Moreover, 
it does not describe correctly the evolution at very short distances. In order to 
improve we perform the second step of our program: the correcting function 
AA^ is introduced to incorporate the correct DGLAP kernel at short distances. 
For AA^ we write down a DGLAP-type linear evolution equation. 

The non-linear evolution equation for the function N has the form: 
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The above equation was derived in quite different approaches to high density 
QCDEI and hence we believe it to be a very reliable tool. The equation (|l|) 
is derived for the large number of colors Nc and for a constant as ■ We have 
dropped the 6-dependence implying the large b limit. Having assumed b to 
be large we will allow to ourselves to extrapolate back to 6 = 0. The non- 
linear equation (]l|) is a subject to certain initial conditions, which we set at 
xo = 10^^. The initial conditions are taken in the Glauber-Mueller form. 



3 Numerical solution 

The equation (|^) is solved numerically by the method of iterations^ for 6 = 0. 
The solutions are obtained for the fixed Us — 0.25 and for the one loop 
running as- The figure |l] presents examples of the solutions as a function 
of distance. The solutions display a step like behavior consistent with all 
theoretical predictions. At short distances the function N tends to zero, while 
at large distances it approaches unity, which is the unitarity bound. 
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Figure 2. Saturation scale as a function of Igx. 
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Figure 3. Gluon density as a function of Igx. 



4 Results 

Basing on the above solutions several physical quantities can be determined. 

• The saturation scale Qs{x) is estimated as a typical scale where the 
function N undergoes the step like transition shown above. Though no exact 
mathematical definition of the saturation scale is known we proposed several 
reasonable definitions which in average produce our predictions depicted in 
the figure ^. 

• The solutions found display the scaling phenomena. Namely, the func- 
tion N is not a function of two independent variables x and r±_ , but a function 
of a single variable r = r±Qs{x): N{x,r±) ~ N{rj_Qs{x)). The scaling holds 
with a few percent accuracy in a broad kinematic domain below x =_10~^. 

• The gluon density was defined according to the Mueller formulafl, which 
relates the gluon density to the dipole elastic scattering amplitude. Examples 
of our predictions are shown in the Fig. ^ At small x the glpon density is 
damped by a factor 2-3 comparing to the DGLAP predictionsQ. 

• The structure function F2 is computed. The experimental data below 
X = 10~^ is reproduced with less than 20% error for < 50 GeV^. It is 
important that the data is reproduced without use of the full DGLAP kernel. 

• The LO BFKL equation, which is the linear part of the equation (|l|), 
is solved numerically by the same method of iterations. The obtained results 
are compared in the Fig. ^ with the corresponding solutions of the non-linear 
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Figure 4. Comparison between the solutions of the LO BFKL equation and the Eq. (h 



equation The solution of the BFKL equation rapidly diverges from the 
solution of the non-linear equation. We conclude that the shadowing effects 
become important before the BFKL dynamics actually takes place. 



5 Summary and plans for the future 

A new method for extrapolating the parton densities to high energies was 
proposed. As the first step in our program we solved (numerically) the non- 
linear evolution equation (|^). From the solutions obtained we estimated the 
saturation scale as a function of x. We found an approximate scaling behavior 
of the solutions. Both xG and F2 are found to be significantly damped at high 
energies. We predicted the LO BFKL dynamics unlikely to be ever seen. We 
work currently on single diffractive dissociation from the non-linear evolution. 
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